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Abstract. We examine the influence of non-ideal plasma-density and non-Gaussian 
transverse laser-intensity profiles in the laser wakefield accelerator analytically and 
numerically. We find that the characteristic amplitude and scale length of longitudinal 
density fluctuations impacts on the final energies achieved by electron bunches. 
Conditions that minimize the role of the longitudinal plasma density fluctuations are 
found. The influence of higher order Laguerre-Gaussian laser pulses is also investigated. 
We find that higher order laser modes typically lead to lower energy gains. Certain 
combinations of higher order modes may, however, lead to higher electron energy gains. 
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1. Introduction 

Laser wakefield accelerators (LWFAs) pQ have the potential to play an important role in 
a wide range of applications, including high-energy physics experiments (by providing 
more compact accelerators [21 El IH El E]), nuclear physics experiments (by providing 
electron beams which can be used for the photonuclear activation of nuclei [3 [8]), and 
medical applications [9] (by providing, for instance, compact x-rays sources [TQl EH])- 
For these applications, the recent developments of LWFAs, which currently deliver low 
divergence (~ 10 mrad), high-energy (~ 1 GeV), high charge (~ 10 pC) electron beams, 
with < 10 % energy spread, in externally- [5j or self-guided regimes [5J, are placing 
strong emphasis on the shot-to-shot stability of the accelerators. It is therefore useful 
to understand the influence of shot-to-shot fluctuations on key outputs of LWFAs. 

Shot-to-shot fluctuations of the plasma density and laser intensity profiles impact 
on the reproducibility of LWFAs. In particular, changes in the laser and plasma profiles 
may affect beam divergence, emittance, pointing stability, energy spread, and maximum 
energy. In this paper we will focus on their influence in energy gain. Previous research 
indicated that fluctuations in experimental plasma density profiles may prevent self- 
injection, thus constraining the applicability of LWFAs [12J. However, plasma density 
ramps [13] , typically present in experiments with gas jets, can be used to control the 
phase-velocity of the plasma wave, and facilitate the self-injection of plasma electrons. 
Fluctuations of the laser intensity profile, which frequently exhibits higher-order modes, 
can affect the quality of the electron beams |14j . Nevertheless, properly tailored higher- 
order laser modes can be used to lower the self-trapping thresholds [15J or increase 
the amplitude of the self-injected beam betatron oscillations [151 E]. Thus, either 
to obtain low shot-to-shot fluctuations on the outputs of LWFAs, or to obtain highly 
optimized acceleration regimes, an adequate control of the plasma density [12], and laser 
intensity [H] is required. 

In order to investigate the influence of non standard plasma density and laser 
intensity profiles we use analytical modeling, and particle-in-cell simulations. Standard 
particle-in-cell algorithms are very computationally expensive. Reduced codes such 
as QuickPIC [T5], or advanced PIC algorithms in boosted frames [2"Ul |2"T] are ideal 
to perform systematic parameter scans of the LWFA. In the boosted frame technique 
no physical approximations are involved. The associated computational gains result 
from bringing closer the disparate scales of LWFAs, i.e. the laser and the plasma 
wavelengths. The analysis of the relevant laser and plasma dynamics, however, is 
not straightforward, since simultaneous events in the boosted frame occur at different 
instants of time in the laboratory frame [22] . The use of reduced models can circumvent 
these difficulties. The Quasi-Static Approximation (QSA), for instance, is widely used 
to model LWFAs in reduced PIC codes [231 El]- It precludes the self-injection physics, 
but the dynamics of the laser, and acceleration of externally injected electron beams can 
be investigated systematically. Typically, the plasma response is determined through a 
Quasi-Static [231 El] plasma field solver, and the laser is advanced in the Ponderomotive 
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Guiding Center Approximation (PGCA) [24]. Simulations under the QSA and PGCA 
require lower spatial resolutions in comparison to standard full PIC codes. In addition, 
larger time steps can also be used [251 EH] . Thus, reduced codes employing the QSA and 
PGCA can be much faster than standard full PIC codes. Here we employ the quasi- 
static, massively parallel, fully relativistic, PIC code QuickPIC which can be more than 
two orders of magnitude faster than standard full PIC codes [19] . 

In this paper, we investigate the influence of non-standard plasma and laser 
configurations in the outputs of LWFAs for a wide range of parameters relevant for 
state-of-the-art LWFAs. In Section |2j a reference simulation using state-of-the-art laser 
and plasma parameters is presented and discussed. Then, in Sec. [3j the influence of 
inhomogeneous longitudinal plasma density perturbations is considered. It is shown 
that the impact of the density fluctuations can be minimized when the plasma-density 
variations are small in comparison to the background density, and when the typical 
wavelengths of the perturbations are much shorter than the acceleration distance. 
Section [4] analyses the influence of non-Gaussian transverse laser intensity profiles in 
the acceleration. Analytical estimates for the energy gain are derived and confirmed 
with QuickPIC simulations. A parameter region where higher order modes can lead to 
higher energy gains was also identified. Finally the conclusions are stated. 

2. Self-guiding of 10 J laser pulses 

In this Section, the propagation of a 10 J laser pulse, with parameters close to Ref. [6] 
is investigated using QuickPIC [T9"] . 

In the baseline simulation conditions that we use in this paper, a laser pulse 
with central wavelength Ao = 0.8 /zm, normalized peak vector potential a = 3.85, 
spot-size Wq = 19 /xm (1/e), and duration r = 55 fs (1/e) is initialized at the 
entrance the pre-formed, transversely uniform plasma. The plasma density rises 
linearly to uq = 5.7 x 10 18 cm -3 for 0.065 cm. Uniform density is then used in the 
following 0.73 cm, falling linearly until the end of the plasma, during 0.12 cm. The 
simulation uses a window which moves at the speed of the light in vacuum (c), with 
250 /im x 250 /im x 80 /im, divided into 512 x 512 x 256 cells, with 4 particles per cell. 
An electron beam with Q = 0.16 pC is at the back of the first plasma wave, in a region 
where the accelerating fields are maximum. The charge of the beam is much lower than 
the beam loading charge [27J [28] , in order to provide a quasi test particle acceleration 
regime. The electron beam density profile is given by: 

n «~» { JI if )">{-£)• « 

where £ = z — ct, z is the propagation distance, t the time, £ the center of the electron 
beam, r = a/x 2 + y 2 is the transverse coordinate, o~ z = 1.2 /im is the length of the 
electron beam and o~ r = 1.2 its width. The electron beam was initialized with a 
relativistic factor 7 = 200, with corresponding longitudinal velocity above the phase 
velocity of the plasma wave. Furthermore, the simulations included radiation damping, 
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which results in a decelerating-like force, due to the betatron motion of the beam 
electrons betatron oscillations [29J . For these parameters, we have found that radiation 
damping is negligible [30] . 

The laser pulse peak ao, spot-size Wq, blowout radius r&, and electron beam peak 
energy, are depicted in Fig. [T] For z < 0.1 cm, strong laser spacial and temporal 
compression (Fig. [T^i) occurrs, where the ao increases from a = 3.85 to ao = 13 [31J. 
Simultaneously, the Wq decreases by the same factor, from Wq = 19 /zm to Wq = 7.5 jum 
(Fig. [T]d). As a result, the ponderomotive force F p oc oq/Wq increases nearly by one order 
of magnitude, leading to complete electron blowout after z = 0.1 cm. This behavior, 
also present in [32J, may be due to the mismatch between the initial laser spot-size and 
blowout radius r^. 

For 0.1 cm < z < 0.8 cm, both the ao and Wq variations are smoother (Fig. [T]) 
indicating that matched conditions established self consistently. In average, ao increases 
from ao = 8 to a = 14 until z = 0.8 cm, rapidly decreasing from then on. The increase 
of ao is accompanied by spot-size oscillations between Wo = 8 /iin and Wo = 16 /xm, 
indicating that the body of the laser pulse is self-guided. Figures and [T]d also 
show that, for 0.1 cm < z < 0.8 cm, both the ao and Wq vary by more than 40% 
in each oscillation. The corresponding variations of r&, however, are lower than 10 %, 
thus guaranteeing a stable acceleration regime (in the blowout regime, the accelerating 
gradients are solely determined by r& [33J ) . A closer investigation of Figs. [IJi and 
reveals that Wq and ao oscillate out of phase, i.e. ao is minimum for maximum Wq, 
and vice-versa, which lead to large variations of the ponderomotive force. However, the 
ponderomotive force at r ~ r& is always similar, thus resulting in lower variations for r&. 

The stable accelerating fields lead to beam peak energies of almost 1 GeV at 
z ~ 0.7 cm (in agreement with |6j), which saturated from then on. The stabilization of 
the final electron beam peak energy was a result of the single beam electrons dynamics 
(not from the phase-space rotation due to the dephasing of the electron beam as in [34^): 
a fraction of the beam electrons are in regions of accelerating fields and gain energy. 
The remaining fraction are in a region of decelerating fields, and loose energy. As a 
result, the energy spectrum broadened with constant peak energy. 

These results establish the baseline conditions to analyze possible perturba- 
tions/deviations from the idealized parameters. 

3. Electron acceleration in inhomogeneous plasma density profiles 

Nondeterministic plasma density fluctuations present in experiments impact on the 
acceleration gradients and final energies of electron bunches. Understanding the role 
of density non-uniformities in the LWFA is thus important to design stable accelerators. 
To this end we estimate analytically the single particle energy gain in a LWFA with 
small density perturbations. The expression for the energy gain of a trapped electron 
in the frame that moves with the laser group velocity v g o, and in a reference plasma 
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Figure 1. Baseline QuickPIC simulation results, (a) Evolution of the normalized 
peak laser pulse vector potential ao associated with an initially transversely Gaussian 
laser pulse, (b) Evolution of the laser pulse spot-size (solid line) and blowout radius 
(dashed line) . The blowout radius corresponds to the location of the maximum plasma 
density for a fixed transverse slice, located roughly at the center of the laser pulse, (c) 
Evolution of the electron beam peak energy. 



frequency u p0 , is given by: 

AE = g P ^M d ^ (2) 

J-r b 1 ~~ P<t>0 

where -E acce i is the accelerating field, £ p = z p — v g0 t = z p — v^t ~ z p (l — (3^) is the 
distance of the externally injected electron to a point moving with the phase velocity of 
the plasma wave v^q = P^c = v g0 , and z p = ct is the electron longitudinal position that 
moves at the speed of light. We consider that the electron, initially placed at £ p = — r&, 
i.e. at the back of the bubble, accelerates in the positive £ p direction until the center of 
the bubble £ c is reached. At the center of the bubble, -E acce i = 0, and the acceleration 
stops. 

In order to determine an expression for -E acce i we assume that the center of the 
bubble £ c travels at the wake phase velocity in the blowout regime = c(l — au p /cjQ), 
where a = 3/2, u is the laser central frequency, u p = w p (£ p ) is the position dependent 
plasma wave frequency, and u /u p ^> 1. Since in the blowout regime the accelerating 
fields are linear, we may write -E aC cei = (m e Up/ 'q)((, p — £ c )/2. In order to derive an 
expression for £ c we assume that £ c = when £ p = — r b . Hence the trajectory of the 
center of the bubble is £ c = / Q * accel (^— v^dt, where t acce i is the acceleration time. Since a 
trapped electron moves at c, £ p = (c — v^t, and hence £ c (£ p ) = [co p (^) 2 / lo^q — l] d£. 
In an ideal scenario where the plasma is uniform (u p (^ p ) is constant) the acceleration 
stops at £ c = 0. In general, if Co> p (£ p ) is not constant £ c may be different than zero. As a 
result, Eq. (g) can be re-written as: 

= ~2 J_ n 2^ df - < 3) 

To further evaluate Eq. ^ we assume small plasma density perturbations such 
that = ujp [l + Su p (^ p )} where 5u p (z p ) <C 1 is a small perturbation. For the sake 

of simplicity we also assume that the plasma density varies sinusoidally, and that the 
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average plasma density (frequency) is n (u p o). Hence, 5u p (^ p ) = 5cos(ki^ p + <fi) where 
5 <C 1 is the amplitude of the sinusoidal plasma density perturbation. The initial phase 
= — 7r/2 — ki^ p0 guarantees that 5u p (— r&) = 0, and that the plasma density rises for 
£ p > — r b if 5 > (£ p o = — r b is the initial electron bunch position). The wavenumber k± 
is the typical wave number of the non- uniformity in the frame that travels with t^ - In 
the laboratory frame coordinates (z p (£ p ),t), k\ = (2/3)k 1 i h u' 2 ./uj 2 . 

For the sake of simplicity, we consider that the upper limit of the integral in Eq. ^ 
is £ c = 0. Retaining the leading order terms on the order of 0{5) leads to: 



AE = m e c 



3 u 2 . f^UpQ 38 k p o Uq ( r^Upo k p0 



(4) 



i~ ~ ~WT~^T 1 — _ -^sin(r fe £;i 

'W^o c 2 2 fci uj z p0 \ c ki 

The first term on the right-hand-side of Eq. (jlj), (AE) uni = m e c 2 (3/8)(r b ujp/c) 2 (jo' 2 /oo 2 , 
corresponds to the maximum energy gain in an uniform plasma. Making the substitution 
r b = 2^/a^c/uj p yields (AE) uni = m e c 2 (3/2)a u} 2 l /uj p l Q which is close to the result from 
Ref. [SDj . The differences in the numerical factors are due to the fact that Ref. [3U] 
considers a constant average E acce i. Instead, Eq. Q uses linear accelerating fields. 
The second term, (AE) sin = —m e c 2 (35ulk p0 /2u 2 k 1 )[r b ijj p o/c—(k p0 /ki)sm(r b ki)} are 
corrections associated with the density non-uniformities. These corrections are due 
equally to the variations on the accelerating gradients, and to the variations of the 
laser group velocity through fluctuations of £ c . According to Eq. Q, the ratio 
(A£) sin /(A£) uni is: 

(A£) sin _ A5[r b k x - sin(r b fci)] 

(AE) uni ~ r 2 kl ■ [b) 

Equation ^ is plotted in Fig. [2] as a function of r b k\. The ratio (AE) sin / (AE) uni 
is directly proportional to 5. Larger final energy fluctuations are hence expected in 
plasmas with larger density fluctuations. In addition, the sign of 5 determines the sign 
of (AE) sin . A density increase (reduction) at the start of the acceleration may then 
lead to larger (lower) energy gains in comparison to uniform plasma densities. This 
is associated with the fact that the acceleration gradients decrease as the acceleration 
progresses. Hence, the density fluctuations at the start of the acceleration will be more 
important than those occurring later. The initial density fluctuations then give the 
dominant contribution for (AE) sin . 

Equation ^ also shows that (A£') S i n depends on the typical wavenumber of the 
density non-uniformities. For r b k\ -C vr the density fluctuations can be neglected because 
their amplitude is arbitrary small in this limit. For r b k\ = tt, (A£ , ) S i n /(A£ , ) un ; = 
45 /n ~ 1.275 is maximum. In the laboratory frame coordinates (z,t) this corresponds 
to a wavenumber k!f h = (37r/2)u 2 /(r b ujQ), and to a wavelength X 1 ^ 13 = 27r/A;J ab = 
(4/3)r b u)l/u 2 . The total acceleration length in the matched blowout regime is L accc i ~ 
(2/3)u 2 r b /u 2 [30]. Thus Xf b ~ 2L accel . In the limit r b ki 3> 7T, or equivalently 
yi&b ^ 2L accc i, (A_E) sin — > because the effects associated with the density non- 
uniformities in the acceleration will average to zero. 
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Equation Q assumes that the density variations are sinusoidal, that the wakefield 
amplitude is constant, and that the acceleration distance does not change, and neglects 
pump-depletion and beam loading effects. In order to verify that Eq. ^ holds even 
when these approximations are not verified, we performed particle-in-cell simulations in 
QuickPIC in conditions relevant for experiments. The baseline simulation parameters 
for the laser and plasma are identical to those of Sec. [2] We first investigate the laser 
dynamics in the experimental density profile shown in Fig. |3^i, which was obtained 
from interferometric measurements of the plasma channel profile during the experiment 
reported in [B]. Figures [3|d, and|3]3 show that the global evolution of a and W in the 
inhomogeneous density profile is identical to that described in Sec. [2] the propagation is 
dominated by strong initial self-focusing that leads to the complete blowout of plasma 
electrons, followed by a self-guided propagation stage. 

In spite of the similar laser dynamics associated with the profiles represented in 
Fig. [3^1, two differences can still be mentioned: first of all, the variation of the laser 
do (Wo) in Sec. [2] is systematically higher (lower) than the corresponding variations 
associated with the experimental profile of Fig. [3^i. Moreover, the frequency of the 
(Wo) oscillations is also higher in Sec. [2] than by using the experimental profile. These 
effects are both due to the higher plasma density used in Sec. [2} On one hand, the use 
of higher densities leads to stronger laser pulse ponderomotive forces, and self-focusing, 
and thus to higher peak a 's. Consequently, the frequencies of oscillation of both W 
and ao are higher. Furthermore, the frequency shifts are stronger and faster for high 
densities, which further increases the peak ao's obtained using the experimental profile 
from Fig. [3} 
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Figure 3. QuickPIC simulation results comparison between an experimental (solid- 
black line) and idealized (dash-black line) plasma density profiles (a) Realistic density 
profile, with small (< 30%) fluctuations from the idealized profile, characterized by 
two ramps at the edges of a uniform plasma density region, (b) comparison between 
the peak normalized vector potentials, and (c) comparison of the evolution of the laser 
pulse spot-size. 



Simulation results from a systematic parameter scan with different amplitudes of 
the plasma density non-uniformities is shown in Fig. |4j The simulated density profiles 
are variations of that of Fig. [3^, with different |5|'s. The value of 5 corresponds to the 
standard deviation of the density variations along the plasma multiplied by a/2 such that 
it is consistent with the amplitude of the sine function used in the analytical model. In 
addition, the average plasma density is identical for each simulated profile, and the test 
electron bunches were injected at the same distance from the laser in each simulation. 

Figure |4)d shows the energy gain of external electron bunches as a function of 
the propagation distance for each of the profiles of Fig. [4^,. Figure |4Jd confirms that 
larger plasma density fluctuations lead to larger energy gain variations in comparison 
to uniform plasma densities. In addition, larger |5|'s lead to lower energy gains since 
the first part of the propagation occurs in regions where the plasma density is lower. 
These findings agree qualitatively with Eq. Q. In order to compare the model and 
simulations quantitatively we plot in the inset of Fig. (4Jd the theoretical value for the 
maximum (AE) S i n /(AE) nn i = 45 /ir, and the maximum energy gain for each simulation 
as a function of 5. The simulation results follow the theoretical results for maximum 
(A£) sin /(A£) uni . 

We also analyzed the role of the plasma density ramps in the acceleration gradients. 
Simulations results with ramp lengths identical to those of Figs. |3^ and but with 
different plasma lengths are illustrated in Fig. ^$p. These simulations indicate that the 
energy gain does not depend on the presence of the plasma ramps. Naturally, larger 
plasma ramps will affect final energy gains. However, our results illustrate that in the 
typical experimental considered conditions, electron energy gains are insensitive to the 
presence of the plasma ramps. 

All simulations showed large energy spreads on the order of 100% at the end of the 
plasma. This is due to the fact that initially the wake is not beam loaded. We use low 
beam charges, such that the acceleration occurs in a test-particle-regime. Consequently 
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the acceleration differed for particles at different longitudinal positions. The results 
from these simulations then did not show a clear influence of 5 on the particle energy 
spread. However, it is expected that plasma density non-uniformities will influence the 
energy spread when the plasma wave is beam loaded. 




Distance [cm] Distance [cm] Distance [cm] 



Figure 4. Energy gain by externally injected electrons propagating in non-uniform 
plasma profiles with different amplitudes of the density variations and plasma lengths, 
(a) Simulated longitudinal plasma density profiles as a function of the propagation 
distance, (b) Corresponding energy gain by externally injected electrons. The inset 
in (b) shows a comparison between the simulated energy gain variation after 0.8 cm 
of propagation with the theoretical estimates, (c) Energy gain by externally injected 
electrons using different plasma lengths with ramps identical to those shown in (a). 



4. Role of higher-order Laguerre-Gaussian laser pulses in the LWFA 

The presence of higher order modes in the transverse laser intensity profile is ubiquitous 
in LWFAs. Hence, understanding the role of higher order modes in the laser propagation 
and in the final electron bunch energy is important to plan and interpret experiments. 
In this section we will examine LWFAs excited by radially symmetric laser pulses with 
higher order Laguerre-Gaussian modes [35]. In vacuum, the slowly varying envelope of 
the normalized vector potential associated with an arbitrary laser beam profile can be 
written as a sum of Laguerre-Gaussian modes as: 

x exp \ i 2 M z \ ) exp (*^) ' ( 6 ) 

where f(z) is the longitudinal laser profile, is the azimuthal coordinate, Wq(z) = 
Wq^JI + z 2 /z 2 is the vacuum laser spot size as a function z, z r = k^W^j2 is the Rayleigh 
length, R(z) — z (1 + z 2 /z 2 ) is the radius of curvature of the laser wavefronts, and l)p 
are generalized Laguerre polynomials of degree p. Since we analyze the propagation of 
radially symmetric laser beams we consider I = 0. Eq. d6]) then becomes: 




Influence of realistic parameters on state-of-the-art LWFA experiments 10 

where L p = L P °K Equation |7) shows how the different combinations of higher order 
Laguerre- Gaussian modes lead to different radial laser intensity, energy and power 
distributions. Since the laser power P is conserved for each transverse laser slice, we 
assume that initial laser power distribution determines the electron energy gain. The 
energy gain AE in the blowout regime as a function P is |30j : 

AE[GeV] = u (nmy> 3 f^!_V /3 (J^Y 3 , ( 8) 

where Ao = 2irc/k Q is the central laser wavelength. 

The fraction of the laser power which is far from the blowout region may diffract, 
not contributing to wake excitation and particle acceleration. Thus, we consider that 
the only laser power useful for particle acceleration is that contained within r < r^. 
According to Fig. [TJthis corresponds to the laser power contained within r < Wo. The 
power at the focal plane contained within a laser spot size for lasers with p = (Gaussian 
beam) and p = 1 is: 



P oc / r 



w r / r 2 \ / 2r 2 \ "I 2 







« (0) + « (1) ex P (- w J L pyw , 



dr 



To make quantitative predictions we assume that the laser pulse energy E User is 
constant for any combination of and According to Eq. ([7]), the laser energy is 
proportional to E la " ser oc (a^ 2 + a^ 2 ) = ag, where ao is the normalized vector potential 
of a purely Gaussian laser beam with Wo and f(z) identical to the higher order mode. 
Thus, a(°) = al — aW 2 . Making use of the latter when inserting Eq. (JoJ) into Eq. (jij), 
we can calculate the ratio AE^ / AE^ given by 

AEC) ~ \ a 2 [e 2 - 1] J ' ( } 

where A£ (1) (AE^ )) is the energy gain of the laser pulse with higher order (purely 
Gaussian) modes. Note that when the higher order laser pulse power inside r < Wo 



is the same as a Gaussian pulse then AE^ / AE^ = 1. Equation (10) can be used 
to examine the impact of the first p = and p — 1 Laguerre-Gaussian modes in 
the acceleration of electrons, and can be readily generalized including the presence of 
modes with p > 1. It shows that the presence of higher order Laguerre-Gaussian modes 
generally leads to lower energy gains in comparison to Gaussian laser pulses. This is 
because the laser power is typically distributed around wider radial regions for lasers 
with aW ^ 0. Minimum energy gains occur when aW = -(a /2) v / 2 + \/2 ~ 0.92a 
for which AE^/AE^ = [(e 2 - 3 - 2^2) /(e 2 - I)] 1 / 3 w 0.625. Interestingly, Eq. Q 
also predicts that the energy gain in the range < < ao/\/2 is higher than that 
of a Gaussian pulse. Maximum energy gain AE^ / AE^ is obtained when aS l > = 
(a /2)v / 2^\72 ~ 0.38a for which AE^/AE^ = [(e 2 + 2^2 - 3)/(e 2 - I)] 1 / 3 w 1.04. 
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In this case more laser power is contained within r < Wo in comparison to a purely 
Gaussian laser leading to AE^/AE^ > 1. 

In order to confirm these findings we performed two sets of QuickPIC simulations 
using lasers with higher order Laguerre- Gaussian modes. The baseline parameters (laser 
energy, spot-size, duration, and plasma density) are similar to those presented in Sec. [2] 
In the first set of simulations we analyzed electron energy gain by lasers with p = 
and p — 1. This permitted to make direct comparisons with Eq. (10). Results are 



summarized in Fig. (J5]). Very good agreement was found between the analytical model 



(cf. Eq. (10)) and simulation results. In particular, simulations confirmed that higher 
order Laguerre- Gaussian laser beams can lead to higher electron energies. Simulations 



also revealed that the theory underestimates the energy gain for (a 



(o) a (i) 



(0,±3.85). 



Simulations in this scenario show that a fraction of the laser initially at r > Wq self- 
focuses and enters the blowout region during the propagation. This mechanism increases 
the total laser power inside r < Wq that can be used in particle acceleration, but is not 



accounted for in Eq. (10) 
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Figure 5. Influence of higher order Laguerre-Gaussian laser beam modes in 
the energy gain in the LWFA. (a) Initial central laser lineouts for each simulated 
scenario, (b) Corresponding energy gain as a function of the propagation distance, (c) 



Comparison between the analytical model Eq. ( 10 1 and the simulation results. 



In a second set of QuickPIC simulations we test the model using "realistic" laser 
pulse profiles. The simulated laser parameters are a numerical fit of measurements of 
the focal spot laser intensity distribution from [6j. The simulation uses = 1.9, 
a« = -0.44, a® = 1.14, a® = 0.05, o< 4 ' = 0.46, = 0.13, = 0.33, a<°) = -0.39, 
and Wo = 30 /im. The plasma density, laser duration and laser energy are similar to 
the baseline simulation parameters of Sec. |2j 

Figure [6] shows the main simulation results. For comparison, Fig. [6^, plots 
the Gaussian laser profile associated with the baseline simulation parameters. The 
corresponding plasma density profile is also illustrated in Fig. ^jp. These figures reveal 
that higher order modes form during the propagation even if the initial laser profile is 
Gaussian. This is because higher order modes appear naturally when Gaussian lasers 
propagate in non-parabolic plasma density profiles such as those associated with the 
blowout regime. 
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Figure [6]>d show simulation results using the fit to the experimental laser profile 
including higher order modes. Since the initial laser profile already contains higher order 
Laguerre- Gaussian modes, the laser filaments at z = 0.3 cm are much more pronounced 
than in Fig. [6^,-b. As the laser propagates it self-focuses leading to the generation of non- 
linear plasma waves in the blowout regime. A sharp electron interface that separates the 
bubble from the surrounding plasma is then created. This thin plasma sheet at r = rj, 
bends the laser pulse wavefronts outwards for r > W ~ r b , and inwards for r < W ~ r b , 
which thus confines the central region of the laser pulse inside the blowout region, and 
prevents the entrance of the side filaments inside it. The blowout region thus acts as 
a spatial filter of the laser pulse. Most of the laser light outside the blowout region is 
then lost via diffraction and pump depletion. This process confirms that only a fraction 
of the initial laser pulse power or energy is used to create the blowout and accelerate 
electrons. As a result, the final electron beam energies are lower in comparison with the 
results from Sec. [2j The electron energy gain as a function of the propagation distance 
is shown in the inset of Fig. [6ji. Despite the theory does not account for the details of 
the wake excitation by lasers with higher order modes, and for the laser evolution, it 
agrees very well with simulations. 

The final energy achieved by LWFAs driven by lasers with higher order modes 
is also connected with the peak laser pulse vector potential a 0L = ^ p a!- p \ To 
illustrate this relation, consider two lasers with identical energy. One of them is purely 
Gaussian with parameters matched to the blowout regime, i.e. W c = 2^a 0G , and 
Pg oc {Wqgo-og) 2 = 4oq G . Equation ^ then gives AEq oc a G- In the second laser, the 
intensity profile is redistributed into well defined filaments sufficiently far from each other 
such that they cannot interact through the plasma. Consider in addition that the central 
laser filament spot-size is matched to the blowout radius, i.e. Wql = < Woo, 

and Pl oc (WoLa L) 2 = 4oq L (this condition may be satisfied after some propagation 
in the plasma). The corresponding energy gain is then AEl oc oql- As a result, 
AEl/AEq oc aoL/a G < 1. This reasoning establishes a direct relation between the 
maximum energy gain in LWFAs with the peak a Q of lasers with higher order modes. 
Note however, that this is only a special case of Eq. (10). Still, for the case of Fig. [6j 
ol = 3.18, a G = 3.85, and then AE l /AEq = 0.82 which is very close to simulation 
results, and to the prediction of the refined analytical model (generalization of Eq. ( 10 ) 
up to modes with p = 7). In general, however, the full driver laser intensity distribution 
needs to be considered to predict the energy gain more acuratelly. 

Similarly to Sec. [3j simulations results in this section showed little correlation 
between the energy spread and the considered laser transverse intensity profiles. This 
is also due to the fact that initially the wake is not beam loaded. It is expected that in 
beam loaded scenarios where the energy spreads are lower, higher order modes will also 
affect the energy spreads. Moreover the wake excitation in the presence of higher order 
modes may contrast significantly with that associated with a Gaussian laser. The wake 
wavelength, and amplitude may then vary using non-Gaussian lasers. This may impact 
in the beam loading charge, in the acceleration distance, and hence in the final energy 
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Figure 6. QuickPIC simulation results comparing a scenario using a purely Gaussian 
laser beam (a)-(b), and a laser with higher order Laguerre- Gaussian modes (c)-(d). 
(a) is a 2D slice at y=0 illustrating the laser vector potential after 0.3 cm. The arrow 
indicates the laser propagation direction. The dashed line indicates the position where 
a transverse line-out was taken in (b). (b) Lineout of the plasma density, and laser 
vector potential after 0.3 cm. The dashed curve is the initial laser profile, (c) is a 2D 
slice at y=0 illustrating the profile of a fit to an experimental laser vector potential 
after 0.3 cm. (d) lineouts of the plasma and laser profiles. The laser profile (dashcd- 
line) is characterized by a central filament, with several filaments in the side. The 
inset in (d) shows the energy gain of externally injected electrons as a function of the 
propagation distance. The theoretical prediction is shown by the dark square, the gray 
curve represents the simulation results with higher order laser modes (H.O.L.), and 
the black curve represents the simulation results with a purely Gaussian Laser (G.L.). 
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spread of the accelerated particles. 
5. Conclusions 

The influence of non-ideal plasma and laser parameters has been examined through 
analytical and numerical modeling. A reference simulation in a self-guided propagation 
regime with state-of-the-art parameters served as a baseline to address the most standard 
deviations to idealized parameters. 

An analytical model was developed to predict the energy gain in non-homogeneous 
plasma density profiles. The analysis showed that the effects of the longitudinal 
density fluctuations are minimized when the amplitude of the fluctuations is low in 
comparison to the background plasma density. In addition, the density fluctuations 
can also be neglected when their characteristic wavelength is much shorter than the 
total acceleration distance. The analytical model can be generalized for arbitrary 
density fluctuations by writing the Fourier series of the plasma density, and performing 
the required integrations. However, based on simulation results using realistic plasma 
profiles, we expect that this will not change our conclusions. The theory also ignored the 
laser pulse dynamics (self-focusing, self-steepening, pump-depletion, self-compression) 
except for the laser pulse group velocity. Still, it captured the relevant physics and is in 
very good agreement with PIC simulations in QuickPIC. 

The acceleration of electrons in the wake of drivers with higher order Laguerre- 
Gaussian modes was also investigated. An analytical model was derived to determine 
the electron energy gain in the wake of a higher order laser pulse. The model assumes 
that only the laser power inside the blowout region can effectively accelerate electrons, 
but neglects the details of wakefield excitation. The theory showed that higher order 
Laguerre-Gaussian modes generally lead to lower electron energy gains. This occurs 
because the laser pulse energy is typically distributed along wider radial regions in 
comparison to Gaussian beams. These results then emphasize that better quality focal 
spots lead to optimized acceleration regimes, also in agreement with previous works on 
self-injection [22] • However, the parameters of higher order Laguerre-Gaussian laser 
beams leading to higher energy gains were identified analytically and confirmed in 
QuickPIC simulations. Simulations using transverse laser profiles with larger deviations 
from Gaussian showed that some laser power initially at r > Wo self-focused, joining 
the blowout region. In these situations the analytical model underestimated electron 
energy gain. Typically, however, the analytical and numerical modeling are in very good 
agreement. 

It is important to remark that higher-order Gaussian beams may lead to optimized 
conditions for wake excitation. For instance, the ponderomotive force which is associated 
with the laser-pulse side filaments modifies the properties (thickness and height) of thin 
electron sheet which surrounds the blowout region [33]. Since the accelerating gradients 
are determined by the details of this thin sheet, the use of higher-order Laguerre- 
Gaussian modes may enhance the accelerating fields of the plasma wave. In addition, 
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the use of higher-order Laguerre- Gaussian laser pulses may also lead to more effective 
self-guiding. For instance lasers with transverse super-Gaussian profiles may be more 
efficiently guided in the blowout. 
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